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ON BOAOTAIL BODIES OF REVOLUTION HAVING MINIMIM WAVE DRAG 
By Keith C. Harder and Conrad Rennemann, Jr. 

SUMMARY 


The problem of determining the shape of slender boattail bodies of 
revolirtlon for mi.nimnm wave drag has been reexamined. It was found that 
mi.ninrum solutions for Ward's slender-body drag equation can exist only 
for the restricted class of bodies for which the rate of change of cross- 
sectional area at the base is zero. In order to eliminate this restric- 
tion, certain higher order terms must be retained in the drag equation 
and isoperlmetric relations. The minimum problem for the isoperimetrlc 
conditions of given length, volume, and base area is treated eis an example. 
According to Ward's drag equation, the resulting body shapes have slightly 
less drag than those detennlned by previous investigators. 


INTRODTOTION 


An approximate expression for the wave drag of slender bodies of 
revolution having zero rate of change of cross-sectional area at the 
base was first given by Von Karmm (ref. 1) . By using this expression, 
together with the calcvilus of variations, several investigators (refs. 1 
to 3) have determined minimum- wave -drag bodies for various isoperimetrlc 
conditions. Later, Ward (ref. Ij-) derived the slender-body approximation 
for the drag of bodies with a nonzero rate of change of cross-sectional 
area at the base. 

In reference 5^ Adams considered several minimum-wave -drag problems 
on the basis of Ward's eqixation. In each case he concluded that the 
minimum-drag body had zero slope at the base. This conclusion implied 
that the minimum shapes for Ward's eqmtion are the same as those for 
Von Karmm's. Recently, Parker (ref. 6) presented a different expression 
for the wave drag of slender bodies and showed that the optimum body 
having given length and base area has a finite slope at the base. Clearly, 
this result is not in agreement with that obtained by Adams. 

In the present paper, the problem of determining minimum-drag boat- 
tail bodies of revolution on the basis of linear theory is reexamined 
with particular en^jhasis on the choice of drag equation, isoperimetrlc 
relations, and method of calculating the body shape. The minimum problem 
for the Isoperimetrlc conditions of given length, volttme, and base area 
is treated as an example. 
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DISCTBSION OF MINIMUM-WAVE-DRAG PROBLEM 


Within the approximations of linear theory, the supersonic flow 
past slender bodies of revolution can he represented hy a distribution 
of sources along the axis of the body. The wave drag D of the body 
may be related to the source distribution (ref. 1 ) by 


to 

plj2 


n z 

f’(x)f*U)loge|x - ||dx d| 
j 


( 1 ) 


provided the source distribution f(x) is zero at the nose sind the beise 
(i.e., f( 0 ) = f(Z) = 0) vhere p is the stream density and U is the 
stream velocity. The coordinate system is shown in the following sketch: 


r 



In the slender -body approximation, the source strength is related to the 
body cross -sectional-area distribution A(x) by 


f(x) = — = 2 jtR(x)R’(x) 
dx 


( 2 ) 


and the restriction that f(Z) = 0 in5)lles that either the body is 
closed (R(Z) = 0) or that the body has zero slope at the base 
(R’(Z) = 0 ). 



NACA TN 514-78 


3 


Several inves'tiga'tors (refs. 1 to 3) bave determined minimum- wave - 
drag bodies for various isoperimetric conditions by applying the calculus 
of variations to equation (l) . However,, as a result of the restriction 
that f(^) = 0, these shapes can be considered optimum only for the 
restricted class of bodies having zero rate of change of cross-sectional 
area at the base. 


Ward's Drag Equation 

An equation which does not have the restriction that f(Z) = 0 was 
proposed by Ward (ref. !(•) on the basis of slender-body theory as 


f'(x)f'(|)loge|x - ||dx d| + 


f'(5)loge(2 - 0^ - f^(^)loge (3) 


M is the Mach number. The source strength is 
again related to the body geometry by equation (2) . 

The problem of determining the source distribution which, minimizes 
the drag given by equation ( 3 ) for given isoperimetric conditions without 
specifying the value of f(Z) at the outset is a variable end-point prob 
lem of the calculuis of variations. In appendix A this problem is consid- 
ered for a general type of isoperimetric condition where it is shown that 
if a mathematical, mi.ninmm exists, it satisfies the condition f(^) = 0. 
The significance of the mathematical solution obtained by the variational 
procedure warrants further consideration since the variational procedure 
assumes the existence of a solution at the outset and, consequently, can 
lead only to necessary conditions for the attainment of an extremum. 

Three mutually exclxisive possibilities must be considered: 

( 1 ) A minimum for Ward's equation exists for the class of bodies 
having all values of the slope at the base and satisfies the condition 
f(l) = 0. 

( 2 ) A minimi mi for Ward's equation exists only for the restricted 
class of bodies for which f(Z) = 0. 

( 3 ) Wo minimum exists for Ward's equation. 
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A single example, not satisfying the condition f(Z) = 0 hut having 
less drag than the shape obtained hy the variational, procedure, is suffi- 
cient to eliminate the first possibility. Perhaps the simplest example 
is the cone which, for given length and base area, has less drag than the 


variational minimum (Von Karman’s ogive) for 3 ^ 0.l64. However, 

a more illuminating example is given by the body 


A(x) = 


A(Z) 


1 + 6 log. 


2* 


+ (Z’ - x)loge/l - 


(0 ^ X ^ Z) ( 4 ) 


where A(Z) is the base area, Z' = Z + e, and e is a parameter related 
to the slope at the base R'(Z) by 


R(Z)loge 

R'(Z) = i_ 

2\l + e logg 


For this body. Ward’s eq,uation (eq. ( 5 )) gives the result that, for 
small 6, 


4itD 

pu2 


A(Z) 


l2 


I + e loge —I 
Z’' 


loSp 



Uogg 


XD ‘ 


[rz- 

v«/ 



+ 


0 



which approaches minus infinity as e approaches zero; that is, as the 
slope and curvat'ure at the base both approach infinity. (The mathematiceil 
symbol 0( ) defines the order of a function.) The body shape for 6 = 0 
is shown in the following sketch: 



NACA TN 3^78 


5 



In order to decide between the second and third possibilities, it 
would be necessary to prove the existence or nonexistence of a miniimmi 
solution for Von Karman's drag equation (eq. (l)). Such a study is beyond 
the scope of the present paper; however, it should be noted that the min- 
imum solutions obtained for Von K 4 nnAn's equation have provided a useful 
guide in the search for low-drag shapes- It is in this same vein that, 
later in the paper, Tnlnimuni problems are considered on the basis of a 
different drag equation. Since existence proofs are not attempted, the 
most that can be claimed is that, if a solution exists, it mxist have a 
certain mathematical form. However, in order to avoid the repetition 
of this qualifying remark in the remainder of the paper, solutions 
obtained by the calcultis of variations are referred to as minimum 
solutions . 

The problem under disciission has been previously considered by Adams 
(ref. 5) ■who correctly determined the necessary conditions for a minimum. 
His interpretation of these conditions was that the optimum boattail body 
has zero rate of change of cross-sectional area at the base. However, 
the proper interpretation is that, if a minimum exists, it exists only 
for the restricted class of bodies having zero rate of change of cross- 
sectional area at the base. 


Parker's Drag Equation 

From the preceding discussion, it is clear that Ward's drag equa- 
tion (eq. (3)) cannot be used to determine minimum-drag boattail bodies. 1 
Parker (ref. 6 ) has shown that application of the calculus of variations 
to the drag equation 


Essentially the same arguments can be used to show that Lighthlll's 
drag equation (ref. "j) , -vdiich was derived for slender shapes with discon- 
tinuities in slope, cannot be used to determine minimum-drag bodies with 
corners . 
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4jtP 

oI]2 


rz-pR(z) rz-pR(z) 

/ / f '(g)f '(x)cosh~^ 

(Z - |)(l - x) - 


pR(l)(x - 1) 


dx d| 


(5) 


which he obtained on the basis of linear theory, yields a ndnimum without 
the restriction that f(i) = 0 for the isoperimetrlc conditions of given 
length and baae area. The body shape so determined has a finite slope at 
the base and less drag than the mathematical mlu lT num for equation ( 3 )- 

Equation ( 5 ) contains some higher order terms which are not included 
in the slender -body approximation to the drag (eq. (3)) since Parker did 
not make the slender -body approximation to the velocity potential in the 
derivation. Apparently, the additional terms are necessary in order to 
obtain minimum-drag shapes withoiit the restriction that t{l) = 0. How- 
ever, it shotild not be inferred that equation ( 5 ) necessarily gives a 
better estimate of the drag of bodies satisfying the assumptions of 
slender-body theory than eqtiatlon (3)* Lighthill (ref. 8) has shown 
that the slender-body equations (eqs. (2) and (3)) are fully as accurate 
as the linearized differential equations of motion for sufficiently smooth 
bodies. Consequently, the slender-body results are theoretically equiv- 
alent to those obtained without making the slender-body approximation. 


Isoperimetrlc Conditions 

The isoperimetrlc conditions most commonly considered have been 
those of fixed length, volume, and base area. In order to carry out the 
mathematical, details of determining the source strength which minimizes 
the drag, the isoperimetric conditions must be directly related to the 
source strength. The slugjlest relations would appear to be those given 
by slender-body theory. However, in order to carry out the analysis on 
the basis of -Parker's equation, certain higher order terms miist be 
retained in the isoperimetric relations. In particular, the limits of 
integration in the isoperimetric relations must be the same as those in 
the drag equation. Furthermore, the analysis can sometimes be simplified 
by including certain additional higher order terms in the integrand of 
the Isoperimetrlc relations. 

The relation between the Isoperimetric conditions and soxrrce strength 
used in the example to be treated in the present paper is obtained by 
approximately satisfying the boundary conditions on a cone passing through 
the nose and base. The linear -theory expression. for A'(x) is 



MCA TN 3^78 


7 



where r is on the body surface. Equation (6) is approximately satisfied 
by evaluating the integral on the cone r = R(Z)— . Then, 



where use has been made of the condition f( 0 ) = 0 . Similarly, the 
volume V is given approximately by 


. n Z(l-6) j 

(z - (Z - 0 ^ - ( 9 ) 


In the derivation of equations (8) and ( 9 ) from equation (7), terms of 
the order of 6%ogg6 have been neglected. The slender -body approxi- 
mation to the isoperimetric relations is obtained by equating 5 to zero 
in equations (8) and (9)* 
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Calculation of Body Shape 

When the calculus , of ■variations is applied to the drag equation and 
isoperimetric relations, the resulting source strength for minimum drag 
contains se'veral constants to be determined from the isoperMetrlc con- 
ditions. The calculation of these constants and the body shape can be 
treated independently of the minimization process. 

Since higher order teims ha'VB been retained in the drag equation 
pnri isoperimetric relations, the question arises as to ■whether similar 
terms should be retained in the body-shape calculation. Theoretically, 
the inclusion of these terms does not affect the accuracy of the result 
for shapes satisfying the assun^itions of slender -body theory. Even so, 
it is interesting to con^tare the various body shapes obtained from the 
source distribution 


f(|) = K 


i\i + 3R(i) - i 


(10) 


found by Parker to give minimum drag for the isoperimetric conditions of 
gi^ven length eind base area. In eqiiatlon (iD) , K is a constant to be 
determined from the isoperimetric conditions. The num er ical value of K 
depends on the method -used ■to calculate the body shape. 

Parker calculated the body shape from this soxrrce distribution -with- 
o\it making the slender -body approximation by numerically solving the 
integral eqxiation 

:tR2(x) = J • f’U)0x - 1)2 - 32 r 2( x ) d| (U) 


The body shape calculated from equation (lO) by means of the slender -body 
expression A*(x) = f(x) is 


A(x) = 


02 


2(1 + c )2 


t|^ 1 - t2 + cos“^(-t) 


(-l^t^c) (12) 



K 
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vhere 


X _ 1 + t 
I 1 + c 


Kl^ ^ a£D 

2(1 + c) c\J± - + cos“^(-c) 


^ ^ I - PR(z) 
I + 3R(l) 


When the expression A'(x) = f(x) is altered to take partial l y into 
account the fact that a given point on the body is influenced only by 
sources in the upstream Mach cone by equating 


dA(x) _ 
dx 



(13) 


the body shape is given by 


A(x) = 


1% 


4c(l + c) 


+ cos'l(-t) 


(-1 ^ t ^ 2c - 1) (l4) 


where 


t = 


2cx 

~T~ 


Z% 

4c(l + c) 


A(l) 


2(2c - l)^c(l - c) + cos"^(l - 2c) 


for p 


The body shapes calculated by means of equations (U), (l2), and (l4) 
nil) 


= 0.2 are compared in figure 1. The differences between the 


I 


shapes are small even for this rather large value of 


For smaller 
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differences are even less and the shapes become 

approaches zero.- Similar body-shape calculations 

and (13) for the source strength that is derived 
in the example (fixed volume, length, and. base area) were performed for 

B = 0.05 and several values of,' where V is the volume. 

I .. 

This comparison is not presented since the body shapes obtsiined by the 
two methods are eOmost identical- Evidently, the difference in body 
shape is appreciable only for shapes that cannot be considered slender. 
Consequently, the sinipler slender -body relation is preferable. 

The discussion concerning the inclusion of higher order terms is 
briefly summarized as follows: Higher order terms must be retained in 

the drag equation in order to obtain the mlnimian-drag boattail body; 
having done this, higher order terms must silso be retained in the iso- 
perimetric relations in order to perform the analysis. Once the source 
strength for tti-^ ninriTm draig has been determined within several undetermined 
constants, higher ordLer terms need not be retained in the calculation 
of the shape and drag of bodies satisfy in g the assungjtions of slender- 
body theory. 


values of p ^ , the 

1 

coincident eis p 

I 

based on equations ( 2 ) 


PROBLEM OF LENGTH, 'VDLLME, AND BASE AREA 


The problem of determining the body shape that gives minimum wave 
drag for fixed length, volume, and base area is treated in order to 
illustrate the ideas developed in the preceding sections. The minimum- 
drag body having given length and base area or given length and volume 
can be obtained as special cases of the problem rinder consideration. 

The source distribution for minimum drag is obtained by applying 
the calculus of variations to equations , ( 8 ), and ( 9 )^ and as shown 
in appendix B, this leads to the source distribution 


fU) = (a + b|) 


Z - I + pR(Z) 


(15) 


where a nnd b are constants to be determined from the Isoperimetric 
conditions . 



NACA TN 5^78 


11 


As discussed in the previous section, the body shape is determined 
on the basis of the slender -body equations. Integration of A'(x) = f(x) 
gives 


( 

CVJ 

II 

A 

t\/l - t2 + cos“^(-t^ - ^1 - 

(1 + c)^ 

2 

L ^ 5 


(-1^ t ^ c) 


(16) 


where 


X _ 1 + t ^ _ I - Q-R(l) ^ _ 1 + mil) ^ 

Z 1 + c Z + PE(Z) 2 


and 


A = 


a + 


I + mil) ^ 
2 


The base axea is given by 


A(Z) = 


Z2 

(1 + c)‘ 


A 

2 


:\/l - c2 + cos"?'(-c) - §(l - c2)^/^ 

^ 5 


(17) 


The voltmie is obtained from equation (16) as 


1 

A 

c cos“l(-c) + \/l - c2 - i(i - 

2(1 + c )5 


V 5 


^1 - c2)^/^ 

5 




C‘=^ + cos 


l(-c) 


(18) 



12 


MCA TN 3^78 


Eq-uation (16) for A(x) (with the use of eqs. (17) and (18)) reduces 
to the mlnimuiii-drag body shape given by Haack (ref- 3 ) Adams (ref. 5) 
when c = 1 . In this case, A(x) is given by 


A(x) 



_ 8 _ 2V 


A(l) 


(1 - t2)5/2 


(_1 ^ t ^ 1) (19) 


where 


X _ 1 + t 
I ~ 2 


This body was obtained by Haack on the basis of Von Karman's drag equa- 
tion (eq. (1)) and by Adams on the basis of Ward's drag equation and is 
referred to in the remainder of the paper as the Haack- Adams body. 

The constants A and B in the equation for the body shape 
(eq. (16}) are determined from eqmtions (17) and (18) in terms of I, 
p, nil), and V. The solution may be expressed as 


= Ai - App2 ^ (20) 

l^ 

and 

P^B = Bi - B2p2 (21) 


r( I ) 

where Aj^, A2, and B2 are fimctions of p ■ Values of Aj^, 

^ R(l) 

^ 2 } ® 1 ^ ®2 given in table I for values of p ^ - between 0.01 

I 

and 0-10. 
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A direct comparison of the drag of the body of the present paper 
vith that of the Haack- Adams body is made on the basis of Ward's drag 
equation. From equation (3) , the drag of the source distribution given 
by equation (15) is 


D 1 Z" 


ait 


a2 


pU^ (1 + c)' 


|cos~l(-c)j + 2c\Jl - c2 cos~l(-c) - (1 - c2) 


4AB(1 - c2) 


+ \Jl - c2 cos~l(-c) 




(2 - 5c^)(l - c2) + 


2c\Jl - c 2 (2c 2 - 1) COS"l(-c) + (cOS“l(-c))^ 
2(A + Bc)^(l - c2)loggj^lt(l + c)^ 


( 22 ) 


The drag of the Haack-Adams body is 


D 2I2 


pU- 


2 


9 ^(11 + 52 X/X _ Mil 




(23) 


The form of equations (20), (21), and (22) indicates that 


pU^Z^ 


is a function of p2 X and p jn figure 2, the drag of the Haack- 

Adams body (eq. (23)) and the drag given by equation (22) are plotted on 
a logarithmic scale for several valties of p help orient the 


R(Z) 
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reader, several body sloapes are shown for p = 1, ^ill. = 0.05* The drag 

~L 

given by eq\xation (22) is somewhat less t ha n that of the Haack-Adams body for 
most values of . For example, for p = 0.05 and p2 — = 0.01, 

which represents a fuselage-type shape, the body given by equation ( 16 ) 
has approximately 7^ percent less drag than the Haack-Adams body. Each 

o V 

drag curve begins at a particular value of p^ — > Oj for a given vadue 

2 ^ 


of 


p smaller values of 


p 2 — give rise to negative body areas. 


Rf 2 ) 

For a given value of P ' ■ ' , the slope at the base of the body is 

2 

positive for small values of p 2 ^ and is negative for large values. 

2^ 

The two drag curves become nearly tangent at intermediate values of 

p2 — for which the body slope at the base is near zero. Actually, the 
2^ 

Haack-Adams body must have less drag for this condition since this body 
gives minimum drag for Ward’s equation for the class of bodies which have 
zero slope at the base. 


The value of 



for miniimim drag is obtained, for a given value 


of 



by equating B to zero in equations ( 17 ) and ( 18 ) . 


This 


procedure gives the optinrum body having a given length and base area. 


In figure 3 the body shape of the present paper is compared with 

the Haack-Adams body for p = 0 .O 5 and p^ — = 0.005 and 0.02. 

^ 2^ 

The bodies are plotted to an expanded vertical scale (expanded 5 times) 
to illustrate the differences which for the most part are small. The most 



NACA OIN 5^78 


15 


significant difference occurs near the base for the larger values of 
p2 ^ where the body given by equation (l6) does not exhibit the reflex 
shape of the afterbody characteristic of the Haack-Adams body. 

The effect of Mach nimber on body shape is illustrated in figure 4 
where the optimum shapes ^vertical scale enlarged 2^ times^ for 

= 0.05 anl — = 0.02 are compared for M = and M = \T^- The 

body shape of the present paper exhibits a small dependence on Mach nxm- 
ber, whereas the Haack-Adams body is independent of Mach number. 


CONCLUSIONS 


The problem of determining the shape of slender boattail bodies of 
revolution for minimimi wave drag has been reexamined and the following 
conclusions are indicated: 

1. Minimum solutions, for Ward's drag equation can exist only for 
the restricted class of bodies for which the rate of change of cross- 
sectional area at the base is zero. 

2. In order to eliminate this restriction^ certain higher order 
terms mvist be retained in the drag equation and Isoperlmetrlc relations. 
However, higher order, terms need not be retained in the calctilatlon of 
drag and body shape from the source distribution. 

5 . Adams in NACA TN 2550 correctly determined the necessary conditions 
for a minimum for Ward's drag equation. His interpretation of these con- 
ditions was that the optimum boattail body has zero rate of change of 
cross-sectional area at the base. However, the proper interpretation is 
that, if a minimum exists, it exists only for the restricted class of 
bodies having zero rate of change of cross-sectional area at the base. 
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4 . Application of the ideas expressed in conclusion 2 to the minimum 
problem of given length, volume, and base area led to body shapes which 
have slightly less drag than the Haack-Adams body. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , June 8, 1955 • 
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APPENDIX A 

APPLICATION OF CALCULIB OF VAEIATIONS TO WARD'S DRAG EQUATION 


In this appendix the calculus of variations is applied to Ward's 
drag equation (eq. (3)) a general type of isoperimetric condition 
to determine the soxirce strength for minimum drag when the source strength 
at the base f(l) is not specified at the outset. From equation (3), 


f '(x)f '(|)loge|x - I dx d| + 


f'(l)loge(^ - 1)11 - f^(^)loge P (Al) 


conditions considered may be related to the 
source strength by expressions of the form 


li = r fU)giU)d| (A2) 

0 

For example, g(|) = 1 for fixed base area and g(|) = (l - |) for 
fixed volume. In the subsequent analysis, it is assumed that one of the 
isoperimetric conditions is that of fixed base area. This assumption 
simplifies the analysis without restricting its generality. 

In. the derivation of eq\iation (Al) it is assumed that 


Ward's drag equation is 



The usual isoperimetric 


f(0) = 0 


(A3) 


Equation (A3) gives one of the end-point conditions to be satisfied by 
the minimizing source distribiition. At the other end point, x = 1, the 
value of f(x) is not prescribed. 
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The source dlstrlhutlon for TnlnlTTrum drag is obtained by considering 
the variation of the function 



(a4) 


where the Lagrange nnjltipliers are determined from the isoperimetric 

conditions. The variation of equation (a 4-) is obtained by considering the ■ 
one-parameter family of con^aris’on functions (see ref. for example) 


f(x) = F(x) + €Ti(x) 


(A5) 


where F(x) is the function which minimizes equation (a 4-) , e is the 
parameter of the family, and ti(x) is an arbitrary function within the 
condition 


T](0) = 0 


(a6) 


This condition arises since all the comparison functions must satisfy 
the same end-point condition as the minimizing function. Since no end- 
point condition is prescribed at x=l, t|(Z) is arbitrary. 

From equations (a 4-) and (A5) , J is a ftinction of e and the source 
strength for minimum drag is determined from the condition 



Tl(^) 


lim 

L 


F(Z)loge 


- i) 

^R(^) 


(AT) 


where the equation has been sin^jlified by severeil integrations by parts 
and use of equations (A5) and (a6) . 
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Since equation (AT) must hold for all choices of ti(|) consistent 
with eqmtion (a 6 ) , it m-ust in particular hold for those choices of ii(|) 
for which ti(1) =0. For such rj(g). 



Tl(l) 


d r ^ F(x)dx ^ 
d| Jq X - I 






0 


and from the basic lemma of the cadculus of variations (ref. 


8 ), 


r ^ F(x)dic 
d| 0 X - I 



(A 8 ) 


With this result, and for general t](|) once again, that is, not 

necessarily eqml to zero, the end-point condition obtained is that 


%fr 


= 0 


(A9) 


In order to satisfy this condition, F(l) must equal zero. Consequently, 
the body shapes which give a ^mathematical minimum for Ward's equation, 
if they exist, must have zero rate of change of cross-sectional area at 
the base. 
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APPiaroix B 

APPLICATION OF CALCULIB OF VARIATIONS 
TO PARKER'S DRAG EQUATION 


The source dlstrlhiitioii for mininrum drag for the Isoperlmetrlc con- 
dition of given length, volume, and, haae area is obtained, by consid.ering 
the variation of the function 


J = 


to 

pl]2 


•t ?V2A(2) 


+ A2V 


(Bl) 


■where D, A(2), and. V are gi-ven by eqtiations (5), (8 ), and. (9)^ respec- 
tively, and. and. A2 are Lagrange multipliers. By proceeding in the 

same manner as in appendix .A, the variation of equation (Bl) is obtained 
by considering the one-parameter family of comparison functions 


f(x) = F(x) + 6Ti(x) , (B2) 

where F(x) is the function which minimizes equation (Bl) , e is the 
parameter of the family, and 'n(x) is an arbi'trary function within the 
condition ti(o) = 0. 

The source strength for miniTmm drag is then determined from the 
condition 


dJ 
dc cbO 


s;. =0 


-X 


i-ps(z) 

0 d| 


2 

'J n 


- if - B^(t) 




° ^ V(i - 


dx + 


’'i + - 5) 


\/(l - t)2 - P%2{1) 





r , . 1 

nU - pB(l)) Um 

, ri-pBU) -.1)^ ^ . L . 

\/(i - D® - 

S— »l-PB(I) 

° ' VG - •- ^ J 



(B3) 
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•where the eq-uation has been sin^ilified by several integrations by parts 
and use of the conditions 
*•* 

F(0) = Ti(0) = 0 

Since equation (B3) must hold for all choices of ti(|) consistent 
•with 'n(O) = 0, it must in particular hold for those choices of ti(|) 

for -which t] |z - = 0. For such 



and from the basic lemma of the calculus of variations. 



where N is a constant. With this result, and for general t]( 0^ "the 
end-point condition is obtained as 
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MCA in 3kj3 


r 


11 m 

I >Z-3R(I) 


<2 


P 


I-PR(I) 


g(x) Vu - 1)^ - P%^(Z) 

"" ' ^ '^(l - x)2 - p2r2(x) 


|Ai + - 5) 


\y(l - 1)^ - p2R2(2)| = 0 


(B5) 


Since equation (B^l-) must hold for all values of and in particular 

for I — >1 - 3R(1), from equation (B5)^ N = 0. Hence, Va - if - 32 r 2 (i) 
can be canceled from each term of equation (rIj-) and the following integral 
equation is obtained for the soixrce strength: 


I 


2-pR(2) 


F(x) 
X - i 


dx 


\j{l - xf - p 2 R 2 (l) 


= -^1 - - I ) ( b 6 ) 


The solution of equation (b 6 ) satisfying the condition F(0) =0 is 


F(|) = f(|) = (a + b|)yi I - I + PR(1) 


where a and b' are constants related to ^ 3 ^ and >s 2 * 


(BT) 
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TABLE I 


COEFFICIENTS OF EQUATIONS ( 20 ) AND ( 2 l) 


p2A = Ai - Aap2 X 

(eq. (20)) 

= Bi - B2?2 X 
(eq. (21)) 

Al 

A2 

Bl 

B2 

0.080572 

0.76164 

0.26804 

18.0291 

.064852 

.66854 

.21958 

18.1466 

.050946 

. 57666 

.17551 

18.2780 

.038807 

.46644 

.15590 

18.4255 

.028390 

.59845 

.10121 

18.5919 

.019652 

.51558 

.071557 

18.7809 

.012554 

.23251 

.046443 

18.9975 

.0070607 

.15669 

.026625 

19.2482 

,0031450 

.088845 

.012091 

19.5439 

,00079066 

.052849 

.0031002 

19.8686 

















Figure l.~ Cong)arison of body shapes calculated from source distribution 

given by equation (lO) by various methods for 3 — (l). = 0.2. 
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